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ABSTRACT 

About seventy percent of intermediate-age star clusters in the Large Mag- 
ellanic Clouds have been confirmed to have broad main sequence, multiple or 
extended turn-offs and dual red giant clumps. The observed result seems against 
the classical idea that such clusters are simple stellar populations. Although many 
models have been used for explaining the results via factors such as prolonged 
star formation history, metallicity spread, differential redenning, selection effect, 
observational uncertainty, stellar rotation, and binary interaction, the reason for 
the special color-magnitude diagrams is still uncertain. We revisit this question 
via the combination of stellar rotation and binary effects. As a result, it shows 
"golf club" color-magnitude diagrams with broad or multiple turn-offs, dual red 
clump, blue stragglers, red stragglers, and extended main sequences. Because 
both binaries and massive rotators are common, our result suggests that most 
color-magnitude diagrams including extended turn-off or multiple turn-offs can 
be explained using simple stellar populations including both binary and stellar 
rotation effects, or composite populations with two components. 

Subject headings: galaxies: star clusters: general — globular clusters: general — 
binaries: general — Stars: rotation 

1. Introduction 

Star clusters are usually assumed to be simple stellar populations (SSPs) that all stars 
formed in a short timescale and with the same metallicity. The color-magnitude diagram 
(CMD) of a star cluster should be similar to an SSP isochrone. However, recent observations 
based on the data that are obtained with the Advanced Camera for Surveys (ACS) onboard 
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the Hubble Space Telescope (HST) showed double, multiple or extended main sequence 
(MS) turn-offs (MSTOs), and dual clump of red giants (RC) for many intermediate-age 
star clusters in the Large Magellanic Cloud (LMC), e.g., NGC1751, NGC 2108, NGC 1846, 
NGC1806, and NGC1783. The shape of CMDs near turn- off looks hke a golf c l ub (hereafter 



golf c l ub" shape). One can r efer to many papers such a s iBertelli et al.l (120031) . iPiotto et al. 



mm. 



Piotto et a. 



fcOOTh. iMilone et all (120091. l2010h. 



Mackev fc Brobv NielsenI J2007h: 



Mackev et al. ( 2008 ): lMilone et al.l (l2009yGoudfrooii et al.l ( 120091 ). iBastian k de MinkI (l2009f l. 
Girardi et al.l (l201ll ). and iRubele et al.l ( 120111 ) for more details. 



Different mode ls ha ve been used for inter preting the special CMDs. For example 
Mackey et al.l (120081 ) and lGoudfrooij et al.l (120091 ) explained this by spread of chemical abun- 



danc e , but it is disagreement with the similarity of na e tallici t y of all stars (iMucciarelli et al. 



2008 



Goudfrooi] et al. 


2011a 


b). Then 


Mackev et al. 


(2008 


); 



preted this by capture of field stars. This reproduces extended MSTO but h as difficulty 



for ex plaining double or multiple MSTOs (MMSTOs), which was found by e.g.. iGlatt et al. 



(120081 ) . Similarly, a picture based on merger of existing star clusters had been brought for- 
ward (IMackey fc Broby NielsenI l2007l ). but it seems not comrn on that clusters in a normal 
molecular cloud merge with each other (iGoudfrooij et al.l 120091 ). A scenario of formation of 
a second generation of stars from the ejecta of firs t generation asympto tic giant branch stars 



was s u ggested but it also has some disadvantages (jD'Ercole et al.ll2008l and lGoudfrooij et al. 



20091 ) . iBastian &: de MinkI (|2009[ ) studied this problem using the effect of rotation on the 
effective and surface gravity of stars, and then on the color and magnitude of star clusters 



his m ode l reproduced t he do uble or multiple population, but the work of iRubele et al. 



(|2010l ) and iGirardi et al.l ( 120111 ) argued that the effect of stellar rotation alone could not ex- 
plain the presence of "golf club" shape of MSTOs. Instead, the y suggested a prolo nged star 
formation history to explain the CMD with "golf club" shape (IGirardi et al.ll201ll ). An age 
spread about 300 Myr is needed for explain "golf club" CMDs including both MMSTOs and 
double RC in such picture. However, this brings a new challenge for the present dynamical 
formation model of star s and star clusters . Some works also studied the possibility of using 
observational selection ( Keller et al.l 2011 ) and un certainty to interpret the observed CMDs 
but it seems not natural (jGoudfrooij et al.ll2011bl ). The kind of model involv ing a mixture 



of sta rs with and with out overshooting ca n not fit the observed CMD, either (jGirardi et al. 



201ll ). A latest model (jPlatais et al.ll2012l ) showed the potential of differential reddening for 
explaining the apparent splitting/widening of MSTOs of an Galactic open cluster, Trumpler 
20, but an age spread is needed. Although most works give their explanation based on single 
stars, a few works a r gue that unresolv ed binary stars may result in extended MS and dual RC 
(IMilone et al.l 120091 : lYang et al.l 1201 ll), but it can not r eproduce some peculia r CMD struc 



tures, e.g., the "golf club" shape (IMackey et al.l l2008l : IGoudfrooij et al.l 120091 : IMilone et al. 
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2009[ ). Even so, the impor t ance of stellar rotation and binaries in modeling CMDs is ob- 



Girardi et al 



2011 



vious flHurlev fc Toutlll998l. iGoudfrooii et al.ll2009l iBastian fc de Minkll2009 



and I Yang et al.ll201l[ ) 



Li et al. 2010 



possibly co ntain a large number of binaries (e.g., Abt 



stars (e.g., iPeterson et al.l 120041 : iMcAlister et al. 



Because it is well k nown that all star clusters 
2006t L and massive rotating 



2005 



9791: 



Lada 



Rover et al.l 120071 ) . it is necessary to 



study the CMDs of intermediate-age star clusters via a combination of the two natural fac- 



Girardi et al. 


2011; 


Yang et al. 


2011) 



previous ones (e.g. 

populations and treatment of rotational ef fect. The effect of stellar rotation is considered us- 
ing the result of Bastian &: de Mink and taking different dis tribution for rotation rate 
cu). The trea t ment of stellar rotation of iBastian &: de MinkI (120091 ) is ac tually different from 
Girardi et al.l ( 120111 ). but it agrees with the results of Geneva group (IBastian fc de Mink 
20091 ). Stars with different masses are considered in our model and a binary fraction similar 
to the observed one is used, which make s the theoretical population more similar to real star 
clusters compared to lYang et al.l (120111 ) . The results are finally found to be different from 
previous works. 

The structure of this paper is as follows. Section 2 introduces the construction of CMDs. 
Then Section 3 shows the main results. Finally, we summarize and discuss on the work in 
Section 4. 



2. Construction of CMDs 



We const r uct synthetic CMDs flowing o ur previous work (e.g., |Li fc Hani [2008al ffl: 



Li et al.l l2010l : IZhongmul l201ll : iLi et al.l |2012|) on modeling binary star stellar population 
(b sSP) and single star stellar population (ssSP). In detail, the initial mass function (IMF) 
of ISalpeterl (Il955l ) is adopted to generate our sample stars. In order to build up binaries, the 
mass of primary component of a binary is gen erated first within the range from 0.1 to 100 
Mq. Note that we do not exclude any stars like lYang et al.l (|201ll ). as two components in an 
interactive binary can transfer and change their masses in their evolution. Then the mass 
of the sec ondary compone nt is calculated via taking a random secondary-to-primary mass 
ratio (g) (IHan et al.lll995l ). q is assumed to obey an uniform distribution within 0-1. The 
separations (a) of two components is given under the assum ption that a is c onstant in logo 
for wide binaries and falls off smoothly at close separation (IHan et al.lll995l ). which can be 
expressed by: 



ania) 



a. 



sep) 



a < ao; 
ao < a < CLi, 



-4- 



where Ogep ~ 0.070, = IORq, ai = 5.75 x 10 i?0 = 0.13pc and m ~ 1.2. This gives about 



half binary stars with orbital periods less than 100 yr (IHan et al.lll995l ) and half single stars. 
After that, a random ecce ntricity (e) within -1 is assigned to each binary, as e affects the 
evolution of stars slightly ( iHurley et al.ll2002[ ). Be cause the binary f raction of real clusters 
is possibly lower than 50% (about 30-40%, see e.g.. lElson et al.lll998l ). we define our sample 
by removing some random binaries from the generated sample. It leads to a new sample 
with binary fraction of 35%, which is then used in this work. 

After the samp le generation, all stars are evolved using the rapid stellar evolution code of 
Hurley et al.l (120021 ) (Hurley code), which uses some formulae fitted from evolutionary tracks 
of stellar models to evolve stars and does not take stellar rotation into account. For stellar 
populations including rotational stars, we add the effect of rotation on effective temperature 
and luminosity to the parameters of massive (> 1.2 M p^) stars that have not lef t main 
sequence (MS), using two fitted correlations presented by Bastian fc de MinkI (120091 ). This 



treatment is followed the work of iBastian fc de MinkI (120091 ). However, we adopt a Gaussian 



distribution with mean and standard deviation of 0.55 and 0. 25 (or 0.15) for rotat ion rate (u), 
according to t he uj distribution of s ome A and F type stars iRoyer et al.l (120071 ) and similar 
application by iBastian fc de MinkI (120091 ) . We also assume that u increases gradually with 
stell ar mass, which gives the fi nal co by multiplying a factor within to 10 for each star (see 
also IBastian fc de MinkI 120091 ) . The factor is linearly increase with mass, and the values for 
1.2 Mq and 1.5 Mq are zero and one respectively. The effect of inclination is rando r nly se t 
to about 0.2 times of that caused by u, according to the result of IBastian fc de MinkI (120091 ). 
A small blue shift of abo ut 0.01 mag in jV — I) of fast rotators is also taken into account, 
according to the result of Platais et al. ( 20121 ). Some limitations in the treatment of stellar 
rotation will be discussed later. Finally, we transform the evolutionary parar neters ([Fe/H], 
T pfj, log q, logL) into colors and magnitudes using the atmosphere library of iLejeune et al. 



(Il998l ). 



3. Results 

Figure 1 shows the CMDs of a few binary star simple stellar populations (bsSSPs) and 
binary star composite stellar populations (bsCSPs) with Z = 0.008. The CMDs in six panels 
are generated by taking different assumptions, and the last four are similar to observed CMDs 
of intermediate-age star clusters in LMC. 

Panel (a) shows the CMD of a bsSSP when all stars could be resolved. It aims to 
show the direct effect of binary evolution or interactions on CMD. A binary fraction of 35% 
and star number of 76 922 are taken for this population. We see that binary interactions 
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lead to some blue stragglers (BSs) and dual RC. BSs are mainly c a.used by mass transfer 



star merger, and possibly high runaway velocity of stars (see also iPols &: MarinusI Il994j ). 
The dual RC results from both normal stars (single and unmerged binary components), and 
merged (most) or interactive binaries. In addition, the turn-off is obviously dominated by 
single stars in the population. This implies that the observed turn-off spread is contributed 
from binary interaction slightly. Note that the turn-off mass in the population is about 1.5 
M0and the deep sequence is close to the isochrone of an ssSSP with the same metallicity and 
age. Furthermore, some red stragglers (or subgiants) (hereafter RSs) are shown on the right 
of MS and under the giant branch. They are primary binary components that transferred 
mass to secondary ones. Some of observed stars in clusters (e.g., NGC 1846) are possibly 
such stars. 

Then panel (6) gives a CMD by taking the spatial resolution of HST ACS into account. 
The distance of LMC is taken as 160 000 light-year, and the angle between line-of-sight and 
the connecting line of two binary components is given randomly. This results in about 85% 
unresolved binaries. We find that MS becomes wider and two parts (blue and red) are 
shown, in which the blue one is obviously dense and narrow. At the same time, double 
turn-offs (2 MSTOs) are presented, but their locations seems different from CMDs of star 
clusters in LMC. Comparing to panel (a), we know that the upper one of two MSTOs is 
caused by unresolved binaries. For convenience, we cite the effect of both binary evolution 
and resolution as binary effect in the paper. Note that the CMD does not have "golf club" 
shape. 

Next, in panel (c) we consider both the effects of binaries and stellar rotation in a 
population. When the effect of rotation is taken into account, the population shows a 
special "golf club" shape similar to the observed CMDs of some LMC clusters. In detail, 
the bsSSP with rotational stars has obvious spread in main sequence and extended turn-off. 
It also shows dual RC and BSs. As a test, we find that the number (about 7) of BSs in 
a fixed range, which is shown by blue lines and affected by field star slightly, is similar to 
that in NGC 1 846 when consideri ng only quality-filtered stars within 30" from the center 



(see figure 3 of lMackey et al.ll2008l for comparison). The distance modulus and color excess 
are taken as 18.45 and 0.1 mag respectively in the comparison. The panel shows that single 
stars, separated binary components, and unresolved stars have some contributions to the 
shape of CMD. 

When comparing the CMD with an isochrone (red points in panel c) of a single star 
simple stellar population (ssSSP) with the same metallicity and age as the bsSSP, we find 
that the ssSSP isochrone locates near the blue and dense MS part of bsSSP. It implies that 
stellar metallicity and age can possibly be estimated by fitting the dense MS using isochrones 
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of ssSSPs. Because the turn-off part lower than the ssSSP isochrone is mai nly caused by 



stellar rotation, it does not indicate more populations. However, as shown by lGirardi et al. 
( 120111 ). it can be explained using prolonged star formation history. 



In panel (d), we show the CMD for the same population, but with a narrow distribution 
of rotation rate for some stars. In this test, half stars are assume to rotate with a Gaussian 
distribution for u, which peaks at 0.55, and with a standard deviation of 0.15. We take 
this assumption because some works show bimodal uj distrib utions that peak near 0.1 and 
0.5-0.6 for some A and F type stars (e.g., iRoyer et al.l 120071 ). Our treatment gives similar 
results as taking two peaks of 0.1 and 0.55 for u. This model reproduces the CMD with 
double turn-offs. Similar CMDs have been observed in some star clusters. Because some 
other works used two ssSSPs to explain similar CMDs, we plot the isochrones of two ssSSPs 
with ages of 1.4 and 1.7 Gyr in this figure as a comparison. We can see that rotation effect 
can also be interpreted as one more stellar population. 

By comparing CMDs in panels (c) and (d) to those of star clusters we find that they 
can fit to the observed CMDs of most intermediate-age star clusters. Besides the extended 
MS or MMSTOs, bsSP with rotation effect can fit the BSs, dual RC, and some RSs (e.g.. 



NGC 1846 and NGC 1987) (see, e.g., iMilone et al.l 120091 ) naturally. Note that RSs can be 
reproduced by assuming more total and resolved binaries. Therefore, it is possible that many 
star clusters are actually bsSSPs with stellar rotation. This supports to the classical picture 
of stellar population of star clusters, i.e., the SSP scenario. 

As some works argued that prolonged st ar formation history may be the best choice for 
interpreting the CMDs of star clusters (e.g.. iGirardi et al.ll201ll ). it is necessary to compare 
the CMD of bsSSP to that of composite stellar populations (CSPs). We do this by panels 
(e) and (/). Panel (e) gives a CMD of a bsCSP with the same metallicity as the bsSSP in 
panel (c). The bsCSP is assumed to form its stars within 300 Myr from 1.5 Gyr by four star 
bursts. We see that this case reproduces wider MSTO and RC. It also has a "golf club" shape. 
Similarly, panel (/) shows a bsCSP with two subpopulations, and considering the effect of 
a narrow u distribution. We find that more MSTOs are generated. This suggests that the 
picture of continuo us star formation history is not the only one for explaining MMSTOs (see 
Girardi et al.ll201ll for comparison). A simple model with two star bursts or a merger of two 
clusters can potentially explain most CMDs with MMSTOs. 
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Conclusion and Discussion 



This paper uses the combined effects of binaries and stellar rotation to interpret the 
CMDs of intermediate-age star clusters in LMC It is shown that most observed CMD features 
can be reproduced via both bsSSP or bsCSP with two subpopulations. In detail, the model 
can show MMSTOs, broad main sequences with two parts, dual RC, and blue and red 
stragglers. It suggests that MS spread and MMSTOs possibly result from both rotation 
and binary effect, while dual RC, blue and red stragglers are mainly caused by binary 
interactions (merger and mass transfer). In addition, the combined effects of binaries and 
stellar rotation leads to CMDs with a "golf club" shape turn-off, which is similar to the 
observed result. Furthermore, bsSPs are shown to have some intrinsic spread in MS, which 
results from unresolved binaries. By taking different assumptions for stellar rotation rate, 
both extended or double MSTO can be reproduced from a bsSSP. Therefore, the special 
CMD shapes of many intermediate-age star clusters in LMC can be explained by including 
both stellar rotation and binary effects in SSP models. In this case, many star clusters with 
special CMDs are possibly SSPs, rather than CSPs. Even if CSPs are needed, some models 
with two components can potentially explain the observed CMDs. 

Although the simple model considering both binary effect and stellar rotation seems 
successful to reproduce the shape of CMDs of intermediate-age star clusters, it is far from well 
understanding the observed CMDs. Firstly, there are many uncertainties in the model. The 
density of blue and red MS parts, the numbers of blue and red stragglers are directly relating 
to the numbers of total and unresolved binaries. Besides, while the treatment of rotation is 
derived from detailed stellar evolutionary models, the effects of rotation on stellar interiors 
and there appearance is considerably uncertain. It should be noted that these models are ID 
models in which rotation is treated in a simple diffusive approximation. On the more massive 
end of t he spectra serious questions have bee n raised about the effects of rotation on interior 



mixing (IHunter et al. 



2011 



Brott et al.ll2011l) . Secondly, the difference between the results of 



this work and iGirardi et al.l ( 1201 ll ) mainly results from various treatment of stellar rotation, 
because different assumptions (single or distributed) for rotation rate can lead to quite 
different CMD features. Although it is reasonable to take a Gaussian distribution for rotation 
rate following the observational result, and assume that massive stars rotate faster since lower 
mass stars have convective envelopes, which may generate magnetic fields through dynamos 
spinning the stars down, our treatment is actually simplistic and artificial. This needs to be 
done better in future studies using either an empirical mass dependence or a grid or detailed 
models accounting for rotation and a certain magnetic braking prescription. In addition, 
two simple fitting formulae were used for modeling the effect of stellar rotation, so the effect 
of rotation on stars' lifetime, which is mainly re sponsible to different color shifts caused 
by rotation in this work and iGirardi et al.l ( 1201 ll ). and the relation between binarity and 
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rotation have not been taken into account. It should be better to use a stellar evolutionary 
code including both binary and rotation in future works. Furthermore, The effects of binary 
interaction and rotation are implemented as two independent effects. In fact, binarity may 
in many cases be the cause of rapid rotation. Moreover, this work did not consider the 
dynamical and chemical evolution of star clusters. They may supply important clues for 
better understanding some gradients for building CMDs, especially star formation history. 
Finally, many factors may affect the observed CMDs, and the roles of them (e.g., prolonged 
star formation and rotation) are somewhat degenerate. This can only be disentangled by a 
series of works involving both the improvement in observation and theoretical modeling. 
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Fig. 1. — Synthetic CMDs. "bsSSP" and "bsCSP" denote binary star simple and composite 
stellar population, respectively. Fraction of unresolved binaries is calculated from binary 
separations and spatial resolution of HST ACS. A^pO is the number of star pairs in a pop- 
ulation at zero age, which is similar for all panels. A^obj is number of objects or points in 
the shown ranges, and for the given populati ons. Two ranges sho wn by blue lines include 
BSs that are compared to those of NGC1846 (IMackey et al.ll2008l ). Errors in V and V — I 
follow Gaussian distributions with standard deviations of 0.01 and 0.014 mag. "nonrotation 
-|- rotation" means half non-rotating stars and half rotational stars. 



